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Abstract. The interaction of a Grand Unification Magnetic Monopole with a nucleon can lead to a barion- 
number violating process in which the nucleon decays into a lepton and one or more mesons (catalysis of 
nucleon decay). 

In this paper we report an experimental study of the effects of a catalysis process in the MACRO detector. 
Using a dedicated analysis we obtain new magnetic monopole (MM) flux upper limits at the level of 
~ 3 • 10~ 16 cm~ 2 s _1 sr _1 for 1.1 • 10~ 4 < \(3\ < 5 • 10~ 3 , based on the search for catalysis events in the 
MACRO data. We also analyze the dependence of the MM flux limit on the catalysis cross section. 
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1 Introduction 

The supermassive magnetic monopoles (MMs) predicted 
by the Grand Unified Theories of the electroweak and 
strong interactions (GUTs) have a core in which the origi- 
nal Grand Unification symmetry is restored |jj . The inter- 
action of a nucleon with the MM core can lead to a reac- 
tion in which the nucleon decays (catalysis of nucleon de- 
cay) . The cross section for this process is of the order of the 
geometrical cross section of the MM core, a ~ 10 -56 cm 2 
||U, practically negligible. In 1982 Rubakov [|||| and 
Callan j(|f7| were independently led to the conclusion that 
this catalysis process may proceed via another mechanism 
and have a much higher cross section, of the same order 
of magnitude of the strong interaction cross sections. 

The Rubakov-Callan mechanism is completely differ- 
ent from the spontaneous decay and so are the correspond- 
ing final states. The main decay channels are e + 7r° and 
/j, + K° for the proton decay and e + ir~ and fi + K~ for the 
neutron decay. According to the theoretical calculations, 
pions should always come with positrons, strange particles 
with muons, and there should be no neutrinos in the final 
states ||. 

The catalysis cross section takes the form [||7],|8| 

where (3 is the MM velocity and ao is a parameter, of the 
order of the cross sections of the strong interactions, which 
depends on the explicit QCD calculations. According to 
some authors |^|| however the cross section for proton 
decay has the form a oc 1//3 2 , so a 1/(3 enhancement fac- 
tor is expected in this case. In addition there could be 
suppression factors for protons in a nucleus || . 

In the hypothesis of a non negligible cross section, the 
MM induced nucleon decay can be exploited to detect the 
passage of a MM in an apparatus. Many proton decay 
experiments (like Soudan pOf , 1MB |ll[] and Kamiokande 
Pjfl ) and neutrino telescopes (such as the Lake Baikal de- 
tector jl3]]) have searched for monopoles via the catalysis 
mechanism: a MM could have been identified as a series of 
decay events more or less along a straight line. All these 
experiments have established limits to the MM flux. 

Beside these direct searches, bounds to the local mono- 
pole flux can be set by means of the catalysis process in an 
indirect way. Astrophysical objects such as neutron stars, 
white dwarfs and also planets can accumulate MMs by 
capturing the ones that impinge upon them [^J. In the 
presence of nuclear matter a MM is able to release energy 



via the nucleon decay catalysis. The limits are established 
by requiring the the energy released by the catalysis can- 
not be greater than the observed one: in this way it is 
possible to set an upper limit to the number of contained 
MMs and then to translate it into a bound to the local 

MM flux [§[n]|i§]. 

MACRO was a large scale multi-purpose detector lo- 
cated in the Gran Sasso underground laboratory, opti- 
mized in the search for GUT magnetic monopoles with 
velocity (3 > 4 • 10 -5 and with a sensitivity well below the 
Parker bound as defined in jL6|. Direct searches for MMs 
were made using various subdetectors both in a stand- 
alone and in a combined way in different ranges of veloc- 
ity. The final results on the upper limits to the MM flux 
are given in Ref. [fl7| ; previous results can be found also 
in [|l8yi9| , |20| , ^l|j22| . These limits refer to a direct detec- 
tion of bare magnetic monopoles of one unit Dirac charge 
(ud = 137/2e), catalysis cross section a ca t < 1 mb and 
isotropic flux. 

In this paper we present a different approach to the 
search for magnetic monopoles using the nucleon decay 
process. We developed a dedicated analysis procedure aim- 
ing to detect nucleon decay events induced by the passage 
of a GUT MM in the MACRO detector; the results of this 
search are reported as a function of the MM velocity and 
the catalysis cross section. After summarizing the prop- 
erties of the MACRO streamer tube subsystem (Sect. 2), 
the subdetector used in the present analysis, we describe 
the simulation technique employed to study the character- 
istics of a catalysis process in the detector (Sect. 3). The 
analysis criteria used in the search for these particular 
events in the streamer tube data are reported in Sect. 4, 
and the results of the real MACRO data analysis in Sect. 
5. The previously published MACRO MM search with the 
horizontal streamer tube subsystem is based on a different 
analysis pattern in which the catalysis process is consid- 
ered negligible; in Sect. 6 we discuss how this analysis 
may be affected by the induced nucleon decay and how 
the obtained upper limits to the MM flux depend on the 
catalysis cross sections. The conclusions are given in Sect. 
7. 



2 The MACRO detector 

MACRO had total dimensions of 76.5 x 12 x 9.3 m 3 , and 
consisted of six "supermodules" arranged in a modular 
structure; each supermodule was divided into a lower and 
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Fig. 1. Ratio of the total number of streamer MM triggers in 
a given sample to the total number of streamer MM triggers 
in the sample without catalysis. 



an upper part, both equipped with streamer tubes and liq- 
uid scintillation counters; in addition one horizontal and 
two vertical planes of nuclear track detector were also 
present. A detailed description of the apparatus can be 
found in Ref.s 23 24|. Special care was taken to ensure 
that the three types of detectors and the readout elec- 
tronics were sensitive to low j3 particles. A single candi- 
date event could have provided distinctive and multiple 
signatures in the three subdetectors, so the experiment 
had enough redundancy of information to attain unequiv- 
ocal and reliable interpretation on the basis of only few or 
even one event. 

We looked for catalysis events in the data collected 
with the streamer tube subsystem. Among the MACRO 
subdetectors the streamer tubes could be the most sen- 
sitive to the catalysis process, having a sufficient combi- 
nation of spatial and temporal resolution to clearly sep- 
arate the MM track from the catalysis products. For the 
same reasons the streamer tube standard analyses could 
be heavily affected by a possible catalysis event, therefore 
we studied as well how the corresponding flux upper limit 
varies with the catalysis cross section. Below we give more 
details on the streamer tube subsystem. 

The lower part of a MACRO supermodule was equipped 
with two horizontal planes of scintillation counters and 
ten horizontal planes of streamer tubes interleaved with 
seven rock absorber layers; the upper part was equipped 
with one horizontal plane of scintillation counters and four 
horizontal planes of streamer tubes without any absorber. 
The vertical west and east walls of the apparatus were cov- 
ered with six planes of vertical streamer tubes interleaved 
with one vertical plane of scintillation counters; the frontal 



walls of the lower part were covered with a similar struc- 
ture. The tubes were filled with a Helium-n-pentane mix- 
ture, in order to exploit the Drell-Penning effect [p5| , p6| 
in detecting MMs, and operated in the limited streamer 
mode. The cell dimensions and the charge drift in the gas 
produced a time jitter, which resulted in a time resolution 
of about 140 ns. 

The streamer tube electronics were designed to be ef- 
fective in the search for slow charged particles with (3 = 
v/c down to 10 -4 . Both the analog and digital part of the 
tube signal were recorded. The analog part was used to 
measure the time and charge of the streamer hits by means 
of a dedicated circuit, the QTP system |^,^4), with a time 
resolution of 150 ns and a charge resolution of ~ 5% the 
typical charge released by a minimum ionizing particle. 
Charge and time information was kept in a cyclic 640 /is 
deep memory. The digital part was used to determine the 
spatial coordinate along a direction perpendicular to the 
tube anode wires, with a resolution of about 1.1 cm; stereo 
pick-up strips were used to determine a second indepen- 
dent spatial coordinate with a resolution of about 1.2 cm. 
Combining these informations it is possible to reconstruct 
the particle track in two space views and in a time view. 

Two independent MM triggers employing the streamer 
tube signals were implemented, one with the horizontal 
planes of the lower part and one with the vertical planes. 
The only assumption made in designing the triggers was 
that a heavy slow-moving particle can cross the detec- 
tor without any apparent variation of direction or speed 
(which is the case of GUT MMs). In this work only the 
trigger on the horizontal planes was used. The primary 
function of this trigger was to measure the time-of-flight 
of particles in the range 3 • 10~ 5 < j3 < 1 and distinguish 
them from random noise. The digital OR of all wires be- 
longing to one plane was summed to the corresponding 
signals from the next module, and then fed to the trigger 
circuitry, which acted as 320 delayed coincidences (called 
(3 -slices) representing the time-of-flight from plane 10 to 
plane 1 with a resolution of 3 fis. The trigger fired when 
the signals from at least seven planes of two consecutive 
modules were aligned in time. 

A complete description of the streamer tube hardware 
and the corresponding triggers can be found in Ref.s (l§|2|||.|7| . 



3 The simulation 

To study the expected characteristics of the catalysis events 
induced in MACRO and their consequences on the data 
analysis and results, a detailed MonteCarlo simulation of 
the physical process in the detector was performed. 

The simulation of the catalysis process was indepen- 
dent of MACRO (and actually can be applied to any de- 
tector set-up). The total cross section for the catalysis 
was left as a free parameter of the simulation; the decay 
channels and the branching ratios were chosen according 
to Ref.s ||]|], namely n,p — > e + ir and n,p — > u + K with 
branching ratios 90% and 10%, respectively. 

Several improvements to GMACRO, the GEANT-ba- 
sed simulation code of the detector, were introduced. They 
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Fig. 2. Ratio of the number of selected events in a given sam- 
ple to the number of selected events in the sample without 
catalysis. 



concern a detailed description of the streamer acquisition 
electronics, which took into account the hit time corre- 
lation when multiple charges on wires and strips were 
recorded, the dead time of the streamer acquisition chain, 
and the charge and time digitization performed by the 
QTP system. It was inserted as well the simulation of the 
streamer tube MM trigger on the horizontal wire planes, 
with a description of its response to particles of different 
velocities and of its interaction with the streamer data ac- 
quisition. In all these tasks we reproduced the behaviour 
of the real electronic components. 

Many samples of simulated events with different (and 
constant) MM velocities and catalysis cross sections were 
simulated. Each sample consisted of 10,000 single MM 
tracks generated from an isotropic flux. Five cross sec- 
tion values were used, namely a cat = 0, 10 -26 , 10 _25 ,5 • 
10~ 25 and 10 -24 cm 2 (cr cat = means no catalysis and 
is used as a reference). For each cross section 5 (7 for 
5 ■ 10~ 25 and 10~ 24 cm 2 ) MM velocities were simulated: 
(3 = 10" 2 , 5 ■ 10" 3 , 10~ 3 , 5 • 10" 4 and 10" 4 (plus 2 ■ 10~ 3 
and 2 • 10~ 4 for the larger a). We kept constant the cross 
section on nucleon a and not the parameter <jq (see Eq. 
|l|) , due to the uncertainties on both the (3 dependence and 
the possible suppression factors on various nuclei. In this 
way our results are independent of the theoretical mod- 
els, and can be used directly to place limits on the cross 
section parameters. 

We did not include in our simulation the rock around 
the detector, unlike other experiments which took into 
account the possibility of a premature trigger due to a 
catalysis event in the surrounding rock. In fact, given its 
trigger logic, the streamer tube MM trigger on the hori- 



zontal planes could not be activated by the decay products 
coming from a catalysis event in the external rock. How- 
ever in case of a fast particle trigger, the trigger signal 
was conveniently delayed to let any slow moving particle 
cross and leave the detector, which was not put immedi- 
ately in dead time [p7j . Therefore we considered only MM 
catalysed nucleon decays inside the detector. 

As mentioned in the Introduction, for the cross sec- 
tion of the induced proton decay there are two theoretical 
models: it could be either proportional to 1//3 as the cross 
section of the neutron decay, or to 1//3 2 . In the latter case 
the proton cross section is enhanced by a factor 1/(3 with 
respect to the neutron decay. To take into account both 
models, for each value of the total cross section a cat and 
for each (3 we performed two simulations: in one of them 
the decay probabilities of a proton and of a neutron were 
equal, in the other one the probability of a proton decay 
was enhanced by a factor 1/(3 with respect to the proba- 
bility of a neutron decay. The results showed a very small 
difference in the MACRO response to the two models, well 
inside the statistical errors. In the following we show the 
results only for the model in which both cross sections are 
equal, the results for the other model being very similar; 
when physical conclusions are drawn, results for both of 
them are given. 

From the experimental point of view it is interesting to 
notice how the detector response changes due to the catal- 
ysis process. Fig. [j] shows the ratio of the total number of 
streamer MM triggers in a given sample over the number 
of streamer triggers in the sample without catalysis (i.e. 
with a cat = 0) for the same (3 (in this Figure and in all the 
subsequent ones the error bars are statistical, and smaller 
than the marker size where they do not appear; only the 
samples at the marker position were simulated and ana- 
lyzed, the connecting lines serve only as a guide to the 
eye). Clearly the number of triggers increases with a cat 
and the MM (3, since the wire signals produced by the 
decay products can contribute to the trigger formation. 
There can be even ~ 40% more triggers for the highest 
values of a cat and (3 which were simulated. The trigger 
increase is small at low velocities and larger at high veloc- 
ities, because the catalysis hits are very near to each other 
in time and can easily contribute to activate the lower trig- 
ger /3-slices (the ones triggered by fast moving particles), 
but can rarely contribute to the higher ones (interested 
by slow moving particles). For a = 10~ 24 cm 2 there are 
so many catalyses per event that the secondary particles 
hits can activate the trigger by themselves, determining 
the larger value of the trigger ratio also for the lowest [3 
value. 



4 Search for nucleon decay events 

All the MonteCarlo samples were analyzed with exactly 
the same programs used later on for the real data. Only 
the horizontal streamer planes of the lower MACRO par- 
ticipated to the trigger formation, but for the event recon- 
struction all horizontal planes belonging to both the lower 
and the upper parts of the detector were used. 
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Fig. 3. Time view of a simulated MM induced nucleon decay inside the detector: the reconstructed slow MM track and the 
fast particles hits (at the arrow position) are clearly distinguishable. 



Events were preselected by requiring 

1. at least one space track matching in both space views 

2. each track must have at least 5 points in each view 

3. in case of multiple wire tracks, at least one "non- 
parallel" track 

4. at least one time track with |ToF| > 0.8 /xs 

Conditions 1 and 2 identify good MM tracks, the third 
one rejects multiple muons from cosmic rays, which some- 
times may confuse the time tracking. The last condition 
selects slow moving particles (a time-of-flight of 0.8 fis 
corresponds to a MM crossing vertically the ten lower hor- 
izontal planes with (3 = 2- 10~ 2 , a value well above the 
upper validity limit of the streamer tube analyses). 

In Fig. the ratio of the number of selected events 
in a given sample over the number of selected events in 
the sample without catalysis is shown as a function of the 
MM (3 and the catalysis cross section a ca t\ this ratio is 
actually the relative efficiency with respect to the case of 
no catalysis. The main reason of inefficiency is a wrong 
ToF estimation because sometimes when there are many 
catalysis hits the time fit may get confused and mix the 
MM signals with the secondary particle ones, giving thus 
a time-of-flight lower than the actual value. Clearly the 
cleaner the event the more correct the ToF computation, 
so for small cross sections this effect is negligible, while 
it is much more relevant for large cross sections. As a 
consequence these criteria result more efficient in selecting 
low a ca t events, and even more efficient for no catalysis 
events. Nevertheless these cuts, the ToF cut in particular, 
cannot be relaxed without contaminating the preselected 
samples from real data with atmospheric muons events 
(especially if accompanied by noise or showers). 

By examining the simulated events it appears that the 
space view of a catalyzed decay is indistinguishable from 
a cosmic muon accompanied by some noise, while it has 
a typical pattern in the time view: a slow track due to 



the MM with some "fast" hits, generated by the cataly- 
sis products, contained in a small ( < 1 fis) time window. 
In Fig. H the time view of a simulated catalysis event is 
shown: the slow track of the MM crossing the detector 
and the hits of the fast decay particles (at the arrow posi- 
tion) are clearly distinguishable . This peculiar signature 
can be exploited in a catalysis-oriented analysis, which is 
applied to the events selected with the preliminary cuts 
described above. The characteristics of the time view of a 
catalysis event, and in particular the presence of fast hits 
grouped together along an almost vertical track, are the 
main feature used in the present search for catalyzed nu- 
cleon decays, while being a possible source of inefficiency 
for the previously published MACRO direct MM searches 
with the streamer tubes, which assumed no decay induced 
by the MM. Although the space view cannot be used to 
tag an event as a catalysis candidate, nevertheless some 
space consistency must be required anyway, for instance 
the catalysis hits cannot be too far from the MM track. 

Since the streamer trigger requires at least 7 planes, a 
track must consist of at least 7 points. In the present anal- 
ysis this condition had to be loosened, since also the catal- 
ysis hits can contribute to the trigger formation, which 
may be activated even by shorter tracks: so each track is 
required to consist of at least 5 points only. If the MM 
traversed the detector, it must be reconstructed also in 
the time view and be a slow particle, hence a (3 cut had 
to be imposed to the time track. There must be a consis- 
tency between the wire and time views: the same planes 
must be present in the wire as in the time fit, since they 
are both reconstructed from the signals coming from the 
same hardware elements; possible hardware inefficiencies 
are taken into account by leaving the possibility of 1 or 2 
planes present in one fit and not in the other. Catalysis 
hits are searched for in ±10 /is around the time track, 
because they cannot be too far in time from the MM. 
At least one catalysis hit must be within 1 m around a 
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Table 1. Results of the real data analysis 
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Fig. 4. Fraction of preselected simulated events with catalysis 
which pass the catalysis specific cuts. 



wire track: the catalysis must have been initiated in the 
vicinity of the MM; in this way we exclude random noise 
hits uncorrelated in space but accidentally near in time. 
All catalysis hits must be inside a 1 (is time window, since 
the decay products are fast particles. Finally, catalysis hits 
must be present on at least 3 out of 4 consecutive planes: 
this is the minimum requirement needed to clearly identify 
a catalysis candidate while rejecting the background. 

The above features have been implemented in the anal- 
ysis of the preselected events following these steps: 

1. select an event with at least one track in both the 
spatial and temporal views with at least 5 points 

2. require no more than 2 planes fitted in either view and 
missing in the other one 

3. for each time track with \(3\ < 10 -2 , which is assumed 
to be the MM track, define a time window equal to the 
time intersection of the track with the first and tenth 
plane ±10 fis in which to search for catalysis hits 

4. if hits on at least 3 out of 4 consecutive planes in a 1 /is 
window are found, of which at least one within 1 m 
from any wire track and the other ones in ±2 m from 
the first one, keep the event as a catalysis candidate. 

Fig. H shows the fraction of preselected simulated events 
with catalysis which meet all the aforementioned condi- 
tions. For the samples at low cross sections {a cat 5s 10 -25 
cm 2 ) this analysis is very inefficient, so our limits will be 
significant only for larger cross sections. The strong de- 
crease of the efficiency for (3 = 10~ 2 is due to the j3 cut 
on the time track. 

The same analysis scheme was applied also to the sim- 
ulated sample without catalysis. No event remained, so 
we can conclude that these cuts have a 100% efficiency 
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in rejecting simulated events without catalysis. As a fur- 
ther check, all events without catalysis contained in the 
samples at lower cross sections were rejected. 

We required that the present catalysis search was valid 
in the range 1.1 ■ 10 -4 < \/3\ < 5 • 10~ 3 , which is the 
same range of the streamer standard analysis: the lower 
limit is the minimum MM velocity required for the Drell- 
Penning effect to take place, the upper limit is required to 
fully reject the atmospheric muon contamination. In Fig. 
[| the fraction of initial events (that is, events contained 
in the MonteCarlo generated samples) which survive to 
all cuts (including the cut to reject the electronic noise, 
described in the next paragraph) for both cross section 
models is shown. This fraction can be considered the total 
efficiency of the present data analysis (the efficiency of 
the hardware is assumed to be the same of the published 
MACRO streamer tube analysis @[l9)). 



5 Data analysis 

The same analysis cuts described in the previous para- 
graph were applied to the whole sample of real data col- 
lected with the streamer tube MM trigger on the horizon- 
tal planes, the same sample used also in [jl?]. These data 
were collected from January 1992 to September 2000, for 
a live time of 71193 hours; the full configuration detector 
acceptance, computed by means of a MonteCarlo simula- 
tion including geometrical and trigger requirements, was 
4250 m 2 sr. 

From a visual scan of a small subsample of preselected 
real events, some of them were found to have heavy elec- 
tronic noise: typically there are two or more horizontal 
planes in which all the QTP channels have fired; for this 
reason the corresponding time hits are biased, so as to 
compromise the ToF measurement. These events repre- 
sent a clearly spurious background which the real samples 
are rich of. In order to reject the noisy events, a further 
cut was imposed on the preselected events before apply- 
ing the catalysis identification analysis, namely to have no 
more than 4 adjacent fired QTP channels. This additional 
requirement had no effect on the simulated samples with 
lower cross sections, and was at the level of few percent 
for higher cross sections. The effect of this cut is already 
included in Fig. ||. 

The events were preselected with the loose criteria de- 
scribed in the beginning of the Section 4; to these events 
the catalysis-specific cuts were applied. In Table 1 the 
number of events which remain after each cut is reported. 
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Fig. 5. Fraction of simulated events with catalysis identified by the analysis with respect to the total number of events in each 
sample as a function of the MM velocity for various cross section values (including the noise cut): equal cross sections for proton 
and neutron decay on the left, enhanced proton decay cross section on the right. 
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Fig. 7. Upper limits to the MM flux as a function of the cataly- 
sis cross section for two different MM velocities, as determined 
by this analysis. 



No candidate survived, therefore upper limits to the mono- 
pole flux were set as a function of the cross section. We lim- 
ited this analysis to the larger values since for a cat < 10 -25 
cm 2 the global efficiency is too low. The time-integrated 
acceptance was the same as for the "standard" streamer 
tube analysis (l7|,|l5| , since exactly the same real data sam- 
ple was analyzed. The new upper limits are plotted as a 
function of /3 in Fig. [| the solid lines represent the stan- 
dard streamer limit in the hypothesis of no catalysis. Only 
the plotted points were actually computed, and the lines 
are only rough interpolations to guide the eye. The lowest 
bounds at /3 = 1(T 4 are 2.99 and 3.38-10" 16 cm" 2 s" 1 sr' 1 
for a cat = 5 • 10~ 25 cm 2 and 10~ 24 cm 2 respectively, 
for the model with equal cross sections, and 2.86 and 
3.36 • 10 -16 cm~ 2 s" 1 sr" 1 for the same cross section val- 
ues, for the model which assumes a proton cross section 
enhancement by a factor 1//3. 

Fig. ^ shows how the upper limit to the MM flux ob- 
tained with this new catalysis-oriented analysis varies as 
a function of the catalysis cross section a cat for two MM 
velocities (3 = 10~ 4 and = 1CT 3 . 



6 Discussion of the results 

As it can be seen from Fig. |g, the MM flux upper limit 
is very good at lower /3, and becomes worse as the MM 
velocity increases. This trend clearly reflects the analysis 
efficiency behaviour. Moreover there quite no difference 
between the two theoretical models on the proton cross 
section form. 

From Fig. |?] we see that the upper limits are very 



stringent for higher cross sections, while for lower values, 
Ocat Si 10 -25 cm 2 they are not very significant. This holds 
for all MM velocities. 

In Table 2 the results obtained by some experiments 
searching for MM induced nucleon decay are summarized, 
together with the MACRO upper limit. It should be noted 
that, except partly for Soudanl, the limits obtained in 
this way rely on the catalysis process, so bounds to the 
MM flux set by these experiments depend on the catal- 
ysis cross section: if the catalysis process does not take 
place there are no limits at all. Moreover the underwater 
and Cherenkov detectors must assume a very high cataly- 
sis cross section to attain the significant upper limits they 
emote. On the other hand MACRO can set significant lim- 
its to the MM flux, at the same level or even better than 
those experiments, with fewer assumptions, and its limits 
stand more or less the same even if the catalysis process 
does not take place or has a negligible cross section. In 
comparison with other experiments MACRO is very sensi- 
tive to small cross sections and can give good upper limits 
even in the presence of higher values of a cat . 

Indirect bounds to the MM flux can be set also by 
measuring the energy emission from astrophysical objects 
like neutron stars and white dwarfs, which could cap- 
ture and accumulate MMs in their interior: these MMs 
could then catalyze the decay of the nuclear matter in- 
side them. The best limit is obtained from the measure of 
the diffuse X-ray background, which essentially the old- 
est neutron stars contribute to, and which is of the or- 
der of KT 23 cm^sr^s- 1 @ (for (3 = ICT 3 and a = 
10~ 25 cm 2 ); limits can be derived from the emission of sin- 
gle massive objects, such as the white dwarf WD 1136-286, 
at the level of ICT 20 cm^sr^s^ 1 fl| (again for /? = KT 3 
and a = 10 -25 cm 2 ). Although very impressive, neverthe- 
less these limits are indirect and very model-dependent, 
especially for what concerns the estimate of the density 
of old neutron stars, based on the assumed birth rate, 
the not well known X-ray absorption in the interstellar 
medium, and the uncertainties in the source distance mea- 
surements. There could be also physical loopholes to the 
validity of these limits [|l4| . Of course if the MMs cannot 
catalyze the nucleon decay there are no bounds at all. 



7 Effects of the catalysis on the standard 
streamer tube analysis 

The standard streamer tube data analysis (along with the 
related hardware) is described in detail in Ref . |1^] . Briefly 
it is based on the search for single tracks in space and time, 
with at least 7 points in the wire and time views and at 
least 6 points in the strip view; all good tracks are required 
to have \(3\ < 5 • 10~ 3 . As in the catalysis search, only 
the horizontal streamer planes of the lower MACRO par- 
ticipated in the trigger formation, both lower and upper 
streamer planes were used for the event reconstruction. 

The standard MM search with the streamer tubes as- 
sumes that the MM does not induce the nucleon decay, 



hence its validity for a cat 5s 10 



" 27 cm 2 



1 mb. To study 
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Table 2. Magnetic monopole flux limits (in cm 2 s 1 sr 1 ) obtained by experiments searching for nucleon decays catalyzed by 
a GUT MM 



Experiment 


Technique 


Flux limit 


(3 range 


SoudanlH 


Proportional tubes 


8.8 • KT 14 


KT* "T" 1 
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1-^3- KT 15 
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Water Cherenkov 
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5 ■ KT 5 ~ 10" 3 
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Fig. 8. Relative efficiency of the standard direct analysis ap- 
plied to the simulated samples with catalysis (using the model 
with equal cross sections) with respect to the simulated sample 
without catalysis. 



the dependency of this analysis on the catalysis cross sec- 
tion, all the simulated samples with catalysis were ana- 
lyzed with exactly the same programs used for the real 
data, applying all but the cut on /3, and the results were 
compared with the simulated sample without catalysis. In 
Fig. U the ratio of the analysis efficiency for the samples 
with a given cross section with respect to the analysis 
efficiency for the sample with a cat = is plotted. The 
behaviour of this ratio as a function of the total cross sec- 
tion is clear: for a cat = 1CP 26 cm 2 = 10 mb the analysis 
efficiency is actually unchanged with respect to the case 
of no catalysis; also for a cat = 1CU 25 cm 2 = 100 mb the 
efficiency is quite the same. So in both cases the catalysis 
does not decrease the MACRO capabilities of detecting 
a MM, and can even help its detection, since there is an 
increase in the number of triggering events (Fig. |l|). In- 
stead for a ca t = 10~ 24 cm 2 = 1000 mb the events are so 
complex and rich of hits that the standard analysis does 
not succeed in reconstructing most of them. This is due 
to the fact that at this level there are so many catalysis 



hits that the reconstruction procedure may fail. We think 
that 1000 mb is the maximum value worth to be consid- 
ered, any simulation for higher values being meaningless. 
From this point of view a cat = 5 • 10 -25 cm 2 — 500 mb 
can be considered an intermediate case: the typical event 
is enough complicated to put in trouble the reconstruction 
code, but not so much as to defeat it completely. 

In addition, the behaviour of the analysis efficiency as 
a function of the MM (3 is interesting: in all samples this 
efficiency is minimum for intermediate velocities, particu- 
larly for (3 — 10~ 3 . At low velocities the time track is very 
inclined, since the MM ToF is large, so the fit procedure 
can more easily distinguish it among the catalysis hits. 
At high velocities the fit may mix the MM hits and the 
catalysis hits, but since the time track is near vertical this 
does not imply a significant error in the ToF reconstruc- 
tion. Instead, at intermediate velocities this hit mixing is 
very deleterious leading often to a wrong ToF estimation. 
In this regime in fact many of the rejected events have a 
reconstructed (3 much greater than the actual one, often 
near to that of a relativistic particle. Clearly, this effect is 
more evident for higher a cat values. 

Since the data set is exactly the same as the one used in 
the standard published analysis, we can use the same ac- 
ceptance and live time, while the analysis efficiency must 
be corrected for the catalysis effects. In Fig. || we show 
the new upper limits to the MM flux as obtained with the 
standard analysis taking into account the catalysis process 
as well (the quoted number is the published limit neglect- 
ing this process). For a cat = 5 • 10 -25 cm 2 the new bounds 



are between 3.43 and 4.83 • 10 -16 cm' 2 s" 1 



sr 



for the 



model assuming the same decay cross section for proton 
and neutron, and between 3.42 and 5.04 • 10 -16 cm _2 s~ 1 
sr^ 1 in the model with enhanced proton decay cross sec- 
tion. 

Fig. ^ shows how the MACRO upper limit to the 
MM flux obtained with the previously published analysis 
changes as a function of the catalysis cross section a ca t for 
two MM velocities (3 = 10" 4 and (3 = 10 -3 . 

From the last Figure it can be clearly seen that the 
MM upper limit from the published standard streamer 
analysis is practically unspoiled by the catalysis process 
for lower cross section values. So the results of the stan- 
dard streamer MM analysis, previously reported for <j ca t < 
1 mb |L7| , can be extended to higher values of the catalysis 
cross section up to a cat ~ 100 mb. Even assuming a cross 
section a ca t — 500 mb, the flux limit would not increase 
more than a factor 2. 
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Fig. 9. Upper limits to the MM flux as a function of the MM velocity for various catalysis cross section, as determined by the 
standard analysis taking into account also the catalysis process, assuming the two theoretical models for the cross section form, 
equal for proton and neutron (left) or proton enhanced (right). 
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Fig. 10. Upper limits to the MM flux as a function of the 
catalysis cross section for two different MM velocities, as de- 
termined by the previously published direct analysis but taking 
into account the catalysis process as well. 

8 Conclusions 

In this paper we presented the results of a dedicated search 
for events of nucleon decays catalyzed by the passage of a 



GUT MM in the MACRO detector. The efficiency of the 
present catalysis-oriented analysis, which used only the 
streamer tube subdetector, is very good for catalysis cross 
sections of the order 100 -j- 1000 mb. Since no candidate 
survived, a 90% CL. upper limit to the MM flux can be 



io- 



sr in the range 



set at the level of ~ 3 ■ 
1CT 4 < (3 < 5 • 10~ 3 . 

We have shown that a dedicated analysis based on 
the nucleon decay catalysis yields flux limits which are 
considerably more stringent than those obtained by nu- 
cleon decay experiments; we also made looser assumptions 
(smaller, more reasonable cross sections, no need to have 
a chain of decay events to clearly identify the MM). We 
were able to set upper limits to the local MM flux even if 
the catalysis process does not take place, using a different 
analysis strategy; here we have shown how this previous 
analysis, developed in the hypothesis of no decay, depends 
on the catalysis cross section. 



Acknowledgements 

We gratefully acknowledge the support of the director and 
of the staff of the Laboratori Nazionali del Gran Sasso and 
the invaluable assistance of the technical staff of the In- 
stitutions participating in the experiment. We thank the 
Istituto Nazionale di Fisica Nucleare (INFN), the U.S. De- 
partment of Energy and the U.S. National Science Foun- 
dation for their generous support of the MACRO exper- 
iment. We thank INFN, ICTP (Trieste), WorldLab and 
NATO for providing fellowships and grants (FAI) for non 



11 



Italian citizens. 



References 

1. Preskill, J. (1984), Ann. Rev. Nucl. Part. Sci. 34, p. 461 

2. Kolb, E.W., Turner, M.S. (1990), The Early Universe, 
Addison- Wesley 

3. Bais, F.A. et al. (1983), Nucl. Phys. B219, p. 189 

4. Rubakov, V.A. (1981), JETP Lett. 33 p. 644 

5. Rubakov, V.A., Serebryakov, M.S. (1983), Nucl. Phys. 
B218 p. 240 

6. Callan, C.G. (1982), Phys. Rev. D26 p. 2058 

7. Callan, C.G. (1983), Nucl. Phys. B212 p. 391 

8. Ellis, J., Nanopoulos, D.V., Olive, K.A. (1982), Phys. Lett 
116B p. 127 

9. Arafune, J., Fukugita, M., Yanagita, S. (1985), Phys. Rev. 
D32 p. 2586 

10. Bartelt, J. et al. (1987), Phys. Rev. D36 p. 1990 

11. Becker- Szendy, R. et al. (1994), Phys. Rev. D49 p. 2169 

12. Kajita, T. et al. (1985), J. Phys. Soc. Jpn. 54 p. 4065 

13. Balkanov, V.A. et al. (1998), Prog. Part. Nucl. Phys. 40 
p. 391 

14. Kolb, E.W., Turner, M.S. (1984), Ap. J. 286 p. 702 

15. Freese, K., Krasteva, E., (1999) Phys. Rev. D59 p. 063007 

16. Turner, M.S., Parker, E.N., Bogdan, T.J. (1982), Phys. 
Rev. D26 p. 1296 

17. Ambrosio, M. et al. Final results of magnetic monop ole 



searches with the MACRO experiment hep-ex 020702C. 

18. Ahlen, S.P. et al. (1994), Phys. Rev. Lett 72 p. 608 

19. Ambrosio, M. et al. (1995), Astrop. Phys. 4 p. 33 and Ref- 
erences therein 

20. Ambrosio, M. et al. (1997), Phys. Lett B 406 p. 249 

21. Ambrosio, M. et al. (1997), Astrop. Phys. 6 p. 113 

22. Ambrosio, M. et al. (2001), [hep-ex/0 110084 accepted by 
Astrop. Phys. 

23. Ahlen, S.R et al. (1993), Nucl. Inst. Meth. A324 p. 337 

24. Ambrosio, M. et al. (1993), Nucl. Inst. Meth. A486 p. 663 

25. Drell, S.D. et al. (1983), Phys. Rev. Lett 50 p. 644 

26. Patera, V. (1989), Phys. Lett A 137 p. 259 

27. De Mitri, I. et al. (1995), Nucl. Inst. Meth. A360 p. 311 



